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Abstract 

A theory is put forward that the electronic phase transition at 0.2 K in Ni- 
doped Bi2Sr2CaCu208 is result of the formation of a spin density wave in the 
system of Ni impurities. The driving force for the transition is the exchange 
interaction between the impurity spins and the spins of the conduction elec- 
trons. This creates a small gap at two of the four nodes of the superconducting 
gap. The effect is to reduce the thermal conductivity by a factor of two, as 
observed. 
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The gap structure of high-temperature superconductors should show up most clearly at 
low temperatures, when the quasiparticles are located near the nodes. Unfortunately, it 
is rather difficult to obtain reliable information in this regime, because only relatively few 
bulk experiments are able to probe properties of the particles and the effect of impurities 
on them. Among these experiments, thermal conductivity k,t offers the most accurate way 
of exploring the DC transport properties. Measurements of kt in Zn-doped YBa2Cus07 
IJhave confirmed the d-wave nature of the order parameter, while experiments in pure 
Bi2Sr2CaCu208 have shown anomalies in a field whose nature has yet to be resolved ||. 

The effect of impurities on kt is expected to be particularly informative because of a 
well-developed theoretical machinery for their calculation ||. It was therefore particularly 
surprising that the doped compound Bi2Sr2Ca(Cui_ x Ni 3 ,) 2 8 has an electronic phase tran- 
sition beginning at a concentration of about x = 0.006 0. As the temperature is lowered 
through T* = 0.2K, kt drops by about a factor of two, indicating a sharp reduction in the 
density of excited quasiparticles. 

The explanation which has been offered for this sharp drop is a transition to a state in 
which the superconducting order parameter breaks time- reversal symmetry, the [d x i_ y i + 
id xy ) state |]5[. This modification is caused by the spin-orbit interaction between the impurity 
spins and the conduction electrons. In this paper I offer an alternative explanation of 
the observations based on the exchange interaction between the impurity spins and the 
conduction electrons which gives rise to a spin density wave (SDW). I shall compare the two 
theories in the conclusion. 

The Hamiltonian is 

n = n d -j^2M(R i )-m(R i ), (i) 

i 

where Tid is the weak-coupling d-wave Hamiltonian and M(R{) is the spin of the Ni atom 
at the impurity site i and fh(Ri) is the spin density of the conduction electrons at the 
same point. This is the usual exchange interaction. It has been taken as local, but this 
assumption is not crucial in what follows. Nor is the sign of J important. The system is 
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treated as two-dimensional. 

I propose that there is spiral magnetic order in the impurity spins below T*: 

(M(Ri)) = M(T) [xcos(Q • Ri) + £sin(Q ■ R { )] . (2) 

Furthermore, take Q = {Q/y/2){l, -1). The reason for this will become clear below. Substi- 
tution of Eq. |^ into Eq. Q shows that there is an effective magnetic field acting on the spins 
of conduction electrons. It contains two terms. There is a coherent field with wavevector Q 
and a random field arising from the fact that the set {Ri} is random. I concentrate on the 
first term here, since it can produce the sudden drop in kt- The effects of the second term 
are discussed below. 

The energy density change in the conduction electron system is 

A£(Q) = - ! l ^ US(k + Q)~ l]Xm(k), (3) 

where S(k) is the static structure factor for the positions of the impurities || and Xm(k) 
is the susceptibility of the conduction electrons. n imp is the two-dimensional density of 
impurities. In order to determine Q, this expression must be minimized. The structure 
factor is that of a highly disordered solid. It will have the peak at k = required by a 
sum rule, and other reciprocal lattice vectors (the first at k ~ y/ni mp ) will be strongly 
suppressed by the Debye- Waller factor. The susceptibility has an unusual structure in a 
d-wave superconductor: Xm(Q) = a Xp£,o\Q\ + • • • + (T / A m )xp(b + c^l\Q\ 2 + . . .) plus terms 
higher-order in T. a,b, and c are model-dependent constants, A m is the maximum gap at 
zero temperature, £o is the coherence length, and xp is the Pauli susceptibility. 

At the very low temperatures of interest here, the non-analytic part is the important one. 
Hence the product AE(Q) has a minimum at \Q \ ~ ^/n impi near the first non-zero peak 
in S(k). Physically, the point is this. The dilute system of impurities can only support 
waves whose wavelength is longer than the interimpurity spacing. The infinite-wavelength 
ferromagnetic state does not lower the energy of the system because the bulk ferromagnetic 
susceptibility vanishes. Because Xm(Q) ~ \Q\, the system chooses the shortest wavelength 



that the dilute impurities will support. The ordering wavevector \Q\ is small, of order 
O.OAA^ 1 . Observe also that if the structure factor is entirely gas-like, then there is no 
transition. There must be some short-range repulsion between the impurities. 

A central point of this paper is that the non-analyticity peculiar to d-wave systems 
is responsible for the transition. This phenomenon is presumably related to the non- 
analyticity obtained in calculations of orbital susceptibilties 0. In s-wave materials, the 
zero-temperature susceptibility x(o) starts with terms of order \q\ 2 : Eq. 1 shows that this 
will strongly suppress the energy of an SDW in a disordered dilute impurity system. Of 
course in a dense sublattice of magnetic atoms the situation is different [|7| and SDW for- 
mation in superconductors has been observed || . 

Considering now the energy AE as an effective interaction in the spirit of RKKY, we 
can calculate the critical temperature of the SDW in mean-field theory for spin 1 on the Ni 
site: 

rp* = J 2 [S(Qo)-i]xm(Qo) U) 
12Ak B fJ.% ' 

Here A is the Taking [S(Q ) - 1] = 0.01, T* = 0.2K, and a susceptibility of fi% x 

l(eV)~ l per conduction electron, we find a very reasonable value of J ~ 50meV for the 
exchange constant. However, because we lack precise knowledge of the values of J and 
[S(Qo) — 1], we cannot use this formula to compare with experiment. Fortunately, most of 
the unknowns in the theory occur only in the product. 

The SDW produces a change in the translation symmetry group. The first Brillouin 
zone is now a thin slice in momentum space. Two of the nodes [ at ±(/cf, kp)/y/2 } are 
in the center of the short axis of the zone and the other two [ at ±(fcp, —kp)/\/2 ] are in 
indeterminate position. Here kp is the length of the Fermi wavevector along the diagonal. 
The states near the former point are strongly affected by the spin ordering because pairs 
connected by the short reciprocal lattice vector have small energy differences. The states 
near the other nodes that connected by the short reciprocal lattice vector have large energy 
differences: this vector is perpendicular to the Fermi surface. 
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To see the effects on the thermal conductivity, we calculate the quasiparticle energies in 
the presence of the coherent part of the effective field in Eq. [I]. The expression for these 
energies near a magnetic zone boundary, obtained by diagonalizing the appropriate 4x4 
matrix, for this is: 

S ±>± (k x , ky) = ± l -{E 2 (k + Q ) + E 2 (k) + 2b 2 eff ± [(E 2 (k + Q ) + E 2 (k)) 2 (5) 
+^b 2 eff (e av (k, Qo) + A 2 av (k, Qo))} 1/2 } 1/2 - (6) 

Here b eff = JxM and E 2 (k) = £ 2 (k) + A 2 (k), where f (fc) are the normal-state quasi- 
particle energies referred to the chemical potential and A(k) is the d-wave gap function, 
while Uik, Qo) = m + Qo) + Z(k)]/2 and A av (k, Q ) = [A(k + Q ) + A(k)}/2. There 
is a level repulsion between the particle-like and hole-like Bogoliubov branches. This 
leads to the development of an energy minigap which is obtained from Eq. |6] by setting 
E(k) = i{k) = A{k) = 0. The result for the minigap energy for the nodes at ±(fcp, kp)/\/2 
is of order bljj-/v s Q , where v s is the slope of the superconducting energy gap at the node. It 
may be estimated as v s ~ A m //cp, where A m is the maximum value of the superconducting 
gap. Then the minigap is roughly lO~ 2 meV when M — 1 . The corresponding gap at the 
other pair of nodes located at ±(k F , —k F )/\/2 is much smaller, of order blff/vpQo, where 
Vp is the Fermi velocity. The big minigap is of order T*, while the latter is smaller by a 
factor of about 30. As a result, the system develops an appreciable gap in two of its nodes 
at the transition, while the other two remain essentially ungapped at T*. 

The thermal conductivity may then be calculated in the presence of M(T), which is 
again obtained from spin one mean field theory. The result from the standard Boltzmann 
equation approach is shown together with the data 0] from a sample with x = 0.015 in Fig. 
1. The curve has been fit from 0.17^ to 0.3K using as adjustable parameters the quantities 
Qo, T*, and b eff with the results Q = 24k B T*/v s ~ k F /50, T* = 0.22K, and b eff = 8k B T*. 
The fact that the fit parameters are in good agreement with the a priori estimates above is 
strong evidence for the correctness of the theory. To the extent that the specific heaat can 
be measured, the entropy change in the same temperature range changes by about iVj mp ln3, 
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which is consistent with this picture. 

The fit in Fig. 1 is poor at low temperatures owing to the neglect of the random term 
in the effective field. It is known that this randomness leads to a crossover from the clean 
result kt ~ T 2 at high temperatures to kt ~ T at low temperatures. This effect is present 
in the data in Fig. 1, though in a temperature regime where the data look noisy. The 
linear regime is not present in the theory as it stands. Improvements in this direction are 
relatively straightforward and will be presented in a longer paper. At higher temperatures 
the data follow a law kt ~ T a with a ~ 1.6 — 1.75, whereas the current theory gives the 
usual result a = 2. There is no theory of this intriguing observation at the present time. 
The corresponding exponent in Zn-doped YBa 2 Cu 3 07 appears to greater than 2 0. 

Considered at the Ginzburg-Landau level, the present theory has the structure: 

F = a M {T)M 2 - a Mm Mm + a m m 2 + 0{M\ M 2 m 2 , m 4 ). (7) 

The order parameters M and m represent the amplitude of the SDW on the impurities 
and on the conduction electrons respectively. The coefficient o>m{T) is entirely of entropic 
origin, implying Om-(T) > 0,as the direct interaction between the impurity spins is negligible 
(/^^imp/^B ~ 10~ 4 K). The temperature dependence of the other two coefficients aum 
and a m may be neglected, as they are proportional to the exchange coupling and the inverse 
zero-temperaure bulk susceptibility, respectively. The transition takes place when the lowest 
eigenvalue of the quadratic form changes sign: this is the equation for T* given above. This 
treatment makes it clear that the SDW involves the ordering of both the impurity spins and 
the conduction electrons. In this regard, the present picture resembles that put forward 
by Balatsky || , which postulates a time-reversal symmetry breaking superconducting order 
parameter (d x 2_ y 2 + id xy ). However, the resemblance is superficial. In the (d x 2_ y 2 + id xy ) 
theory, the Ginzburg-Landau expansion has the 

F = a M (T)M 2 + a \A \ 2 + ai |Ai| 2 + ib MA M z (A At - A AJ) + G(M 2 4 , A 4 ,, A 4 ) 

Here Ao and Ai are the d x 2_ y 2 and d xy superconducting order parameters and M z is the am- 
plitude of the uniform magnetization on the impurity sites. The most important difference 
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from the free energy in Eq. [7] lies in the fact that the coupling term is trilinear. As a result, 
the phase transition can only arise from a change of sign in o,m(T). In view of the estimate 
of the direct magnetic coupling given above, this is unlikely to occur at the temperatures 
in question. The two theories differ also in their predictions for kt- In the SDW theory, 
only one half of the nodes are substantially gapped, leading to a natural explanation of the 
sudden drop by a factor of two. In the (d x 2_ y 2 + id xy ) theory all nodes are gapped. The 
theory must rely on the formation of an impurity band of states in the gap t explain the low 
temperature limit. This would imply that the zero-temperature limit of kt/T is extrinsic, 
depending on the impurity concentration. The fact that kt/T approaches a value compa- 
rable to the universal limit makes this explanation unattractive, but it is clearly consistent 
with the SDW theory. 

In closing it is appropriate to consider how to sharpen the comparison of experiment and 
theory. The two-dimensional SDW transition of course needs to be stabilized by interlayer 
coupling. It is not clear what the threee-dimensional ordering will be, as the dipole or other 
coupling may then compete with RKKY. In any case, the small impurity concentrations 
may make the determination of the structure by neutron scattering difficult. A further 
complication is the dependence of T* on the structure factor. This can cause the transition 
temperature to depend on the preparation method and annealing time. It could be respon- 
sible for the fact that the transition temperatures in samples of Bi 2 Sr 2 Ca(Cui_ x Ni :r .)208 
with nominal compositions x = 0.006 and x = 0.015 are about the same. It is likely that a 
different mechanism suppresses the transition at x = 0.024. At these higher concentrations 
Xmil = 0, T = 0) will become appreciable due to the random interaction term. This will 
favor a transition into a ferromagnetic state, as there is always a 5 2 (k) term present in S(k). 
Repeating the calculation of the quasiparticle energies as above shows that this does not 
produce any gap in the quasiparticle spectrum and thus would not be seen in kt- By the 
same token, application of a uniform magnetic field produces a first-order transition to a 
ferromagnetic state and erases the kt anomaly, as is seen in the experiments. A simple way 
to distinguish between the SDW and (d x 2_ y 2 + id xy ) theories experimentally is that the SDW 
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responds isotropically to an applied field, while the (d x 2_ y 2 + id xy ) would predict the disap- 
pearance of the kt anomaly only if the field is in the basal plane. Furthermore, the SDW 
would lead to transport anisotropy between the (1,1) direction and the (1,-1) directions, if 
domain effects can be eliminated. This could happen automatically if there is coupling of 
the SDW to the superlattice distortion in the crystal structure. 
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FIGURES 

FIG. 1. Comparison of theory(solid line) and experimental data (points)for the thermal 
ductivity divided by temperature versus temperature. The data are taken from Ref. 4. 
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